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Flat-Bands on Partial Line Graphs 
— Systematic Method for Generating Flat-Band Lattice Structures — 

Shin Miyahara*, Kenn Kubo, Hiroshi Ono, Yoshihiro Shimomura and Nobuo Furukawa 
Department of Physics and Mathematics, Aoyama Gakuin University, Sagamihara 229-8558, Japan 

We introduce a systematic method for constructing a class of lattice structures that we 
call "partial line graphs" . In tight-binding models on partial line graphs, energy bands with 
fiat energy dispersions emerge. This method can be applied to two- and three-dimensional 
systems. We show examples of partial line graphs of square and cubic lattices. The method 
is useful in providing a guideline for synthesizing materials with flat energy bands, since the 
tight-binding models on the partial line graphs provide us a large room for modification, 
maintaining the flat energy dispersions. 
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It has been known that several novel phenomena occur in electronic systems with "flat- 
bands", i.e., energy bands with flat dispersions in an entire /s-space. Flat-band ferromagnetism 
is the most well-known example, and many research studies have been devoted to this subject. 
Mielke first proved the existence of ferromagnetism in a Hubbard model on "line graphs". 1-3 
Later, Tasaki showed rigorously the existence of ferromagnetic ground states in the so-called 
Tasaki model. 4 These tight-binding models have a flat-band as the lowest band. Since the 
single-electron states in the lowest band are completely degenerate, the ground states of 
noninteracting models with the half-filled lowest band are highly degenerate due to possible 
spin configurations. The introduction of a Hubbard interaction does not change the energy 
of perfectly ferromagnetic ground states since these states do not encounter on-site repulsive 
interactions. If the energies of other states are increased by such interactions, the so-called flat- 
band ferromagnetism is realized. Such a flat-band ferromagnetism can be stable against small 
perturbations to hopping terms. It is known that a flat-band also appears in a tight-binding 
model on a bipartite lattice with different numbers of sub-lattice sites. The ground state of 
the half-filled Hubbard model on this type of lattice was shown to be ferrimagnetic (Lieb 
theorem). 8 Motivated by these studies, several attempts at finding a flat-band ferromagnetism 
in real materials have been carried out. Following the guiding concept of constructing localized 
wave functions, tight-binding models with flat-bands were proposed in several theoretical 
studies. 9-12 For example, in the graphite with zigzag edges, the edge states formed a flat-band 
and the magnetism caused by the edge band was predicted. 10 
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Imada and Kohno argued that a flat dispersion may enhance the pairing instabilities in the 
systems near a Mott insulator with a singlet ground state. According to them, the enhanced 
degeneracy and suppression of single-particle processes due to the flat dispersion may enhance 
pairing instability in doped systems, which may be the cause of the high-T c in cuprates. 13 

In this manner, the search for tight-binding models with flat-bands is strongly desired 
from both theoretical and experimental points of view. Thus far, several attempts at con- 
structing tight-binding models with a flat-band have been carried out. In most cases, we have 
to determine the parameters in the models as realize a localized wave function by solving the 
Schrodinger equation directly. 11 ' 12 Although a few systematic methods for constructing hat- 
band tight-binding models, e.g., "line graph" and "cell construction", have been proposed, 1, 2,4 
the application of these methods can be possible in certain artificial models. Thus, the number 
of constructed models is limited. In this paper, we introduce a new simple systematic method 
for constructing tight-binding models on a class of lattice structures that we call "partial line 
graphs" in which flat-bands emerge. We can apply this method to any two-dimensional (2D) 
or three-dimensional (3D) lattices. In addition, the generalization of the "partial line graphs" 
is also possible, which makes the flat-band stable in a wide parameter range. Such a flexibil- 
ity is important as a guideline for constructing a material with a flat-band. First, we show 
the method for constructing a partial line graph. Several tight-binding models produced by 
generalizing the partial line graph also realize flat-bands. We describe the rules for the gen- 
eralization and show examples of generalized partial line graphs of square and cubic lattices. 
Finally, we prove the realization of flat-bands in these systems. 

Given a 2D or 3D lattice, we can construct its partial line graphs in the following manner. 
We divide the lattice into N sublattices (N > 1), where each site on a sublattice is equivalent 
and connected only to the sites on different sublattices. Next, we select a site of one (A) of 
the sublattices, which is connected to z sites on different sublattices, and replace the selected 
site with z sites (Ai, A2, • • • , A 2 ) created on bonds connecting the selected site to other sites. 
We connect Ai, A2, • ■ ■ and A z to each other through bonds. Repeating the same process 
on all the sites on the A sublattice, we construct a partial line graph of the given lattice, 
where the A sublattice is replaced by a lattice of clusters composed of z sites. Namely, we 
apply the concept of the line graph only to one of the sublattices. The constructed lattice has 
N + z — 1 sites within a unit cell. A tight-binding single-particle model with the same hopping 
amplitude t on all bonds realizes flat-bands at the energy — t, as will be shown later. 

As an example, let us consider a partial line graph of the square lattice. First, we divide 
the lattice into two sublattices (A and B), as shown in Fig. 1(a). We then create four sites 
(Ai, A2, A3 and A4) on the bonds connected to the sites on the A sublattice (see Fig. 1(b)). 
Connecting the new sites, we obtain a partial line graph shown in Fig. 1(c), which contains five 
unidentical sites in the unit cell. Assuming the hopping amplitude t on all bonds, we calculate 
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Fig. 1. (a) 2D square lattice (z — 4) with two sublattices, where A and B sublattices are shown by 
black and white circles, respectively, (b) All sites on the A sublattice are removed, and new sites 
on the bonds connected to the removed sites are created, (c) The new sites on the clusters are 
connected. The partial line graph of the 2D square lattice is obtained. The unit cell is shown by 
a dashed line. 

the energy dispersion. The result is shown in Fig. 2(a), where doubly degenerate flat-bands 
exist at E = —t. Furthermore, all the dispersion relations along the highly symmetric lines in 
the k-space, k = (0,k y ) or (k x ,0), are flat, as shown in Fig. 2(a). 

As mentioned above, we have assumed that all the hopping amplitudes and the on-site 
potential energies on all sites are the same. However, a flat-band emerges in more generalized 
models containing different hopping amplitudes and/or on-site energies, if we choose these 
parameters properly. The generalized partial line graphs are constructed by the following pro- 
cedures. In the following discussions, we call all the sites other than those on the A sublattice 
simply as B sites, although they may belong to different sublattices. 

(i) The on-site energies on the B sites do not affect flat-bands at E = —t. For example, 
the value of the on-site energies of the B sites e on the partial line graph of the square lattice 
may be arbitrary. We show the energy dispersions for e = 0.5i in Fig. 2(a) together with those 
for e = 0. 

(ii) We can take arbitrary hopping amplitudes on the bonds that connect the A and B sites 
or two B sites, while fixing the amplitudes within A clusters to be t. Flat-bands at E = — t 
are also realized in this case. One of these models is shown in Fig. 3(a) and an example of its 
band dispersions is shown in Fig. 2(b). 

(iii) The hopping amplitudes within A clusters can also be tuned if z > N + 2. We 
classify z sites in an A cluster into M subsets (Si, S2, Sa/), where M should satisfy 
z — N > M > 1. When the hopping amplitudes between sites in the same subset and between 
those in the subsets Si and Sj are considered to be to and %, respectively, the model has a 
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Fig. 2. (a) Dispersion relations for the partial line graph of the square lattice (Fig. 1(c)). The on-site 
energy e on the B site are considered to be and 0.5t, while the on-site energies within four-site 
clusters are 0. (b) Dispersion relations for the generalized partial line graph shown in Fig. 3(a) 
with t\/t = 0.2 and t-ijt = 0.8. (c) Dispersion relations for the generalized partial line graph shown 
in Fig. 3(b). (d) Dispersion relations for the generalized partial line graph shown in Fig. 3(c). 



flat-band at E = —to- For example, we divide an A cluster of the partial line graph of the 
2D square lattice into two subsets, S\ and S2, as shown in Fig. 3(b) and assume to = and 
t\2 = t. The band dispersions of this model with a flat-band at E = are shown in Fig. 2(c). 

(iv) We can even maintain the original sites on the A sublattice, which was removed in 
the original partial line graphs. This type of lattice contains A clusters composed of z + 1 
sites. An example of this generalization of the 2D square lattice and its band dispersions are 
shown in Figs. 3(c) and 2(d), respectively. Note that applying procedure (hi) for the model 
in Fig. 3(c) reproduces the Lieb model. 8 

Thus far, we have shown only the examples constructed from the 2D square lattice. How- 
ever, many types of tight-binding model with flat-bands can be constructed systematically 
by applying the above-mentioned procedures and their combinations to any type of lattice 
structure. As an example of 3D models, we introduce a generalized partial line graph of a 3D 
cubic lattice (Fig. 4(a)) and show its dispersion relations in Fig. 4(b). In this case, flat-bands 
appear at E = and — 2t. 

Let us now prove the existence of flat-bands. We prove first that flat-bands emerge at 
E = — t on an original partial line graph. Electronic states are described using the Bloch wave 
function </>k with the wave vector k: 

0k = Y.R [h k^B 1 (R) H h &JV-1 k0B AT-1 (R) 
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Fig. 3. Generalized partial line graphs of the 2D square lattice, (a) Arbitrary hopping amplitudes 
from the 3i original sublattices. (b) Generalization for the hopping amplitude within the A cluster, 
(c) Addition of the original site on the A cluster. 




Fig. 4. (a) Generalized partial line graph of a 3D cubic lattice, where all hopping amplitudes on the 
bonds are assumed to be t. (b) Dispersion relation. 



+a lk 0Ai(R) + • • • + Wa*(R)] e lk ' R , (1) 

where the wave functions ^ai,-" >^Az represent the atomic states within A clusters, and 
those on B sites are given by 4>bi,-'~ j^BJV-l- The wave function </>k for the energy 
satisfies the eigenvalue equation Ti^k = -E'k'Ak- The Hamiltonian H for the wavevector k is 
represented by an (iV + z — 1) x (N + z — 1) matrix, which can be written as 

nJ B{k) e(k) V (2) 

V ff'(k) T J 

Here, £>(k) and <7(k) are the (N — 1) x (N — 1) and (N — 1) X z block matrices, respectively. 
The zxz block matrix T is independent of the wave vector k, and all matrix elements, except 
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Assume that the eigenfunction <pf^ for the eigenvalue —t has no amplitude on B sites, i.e., 
= for 1 < i < N — lin eq.(l). Then, all z equations (T + tl) t (an c , fl2k> ■ ■ ■ ,i z k) = 
are identical, and the eigenvalue equation is reduced to N independent linear equations for 
the unknown z coefficients aik) a 2k> • ' ' > a zk : 

<?(k) t (aik,a2k,-'- ,a z k) = 0, (4) 

aik + «2kH N 2 k = 0. (5) 

Hence, we obtain z — N independent solutions for z unknowns if the condition z > N is 
satisfied. That is, we have z — iV-fold degenerate energy bands at = —t, whose dispersions 
are flat in the entire k space. For any wave vectors k, by making nodes on B sites, a wave 
function takes the eigenvalue —t. 

It is obvious that the flat-bands at E = — t do not depend on the on-site energy on the 
B sites and the hopping amplitudes from/to the B sites, since they modify only eq. (4) and 
do not change the number of independent equations. Thus, once all the on-site energies in 
the A clusters are identical, we can choose an on-site energy on the B sites and an arbitrary 
amplitude from/to the B sites to realize flat-bands, as was already illustrated for the 2D 
square lattice in Figs. 2(a) and 2(b). 

In generalization (iii), we tune the hopping amplitudes within an A cluster by dividing 
the cluster into M subsets {Si, S2, • • • , Sm} and taking the hopping amplitudes in the same 
subsets to be to and those between the subset Si and Sj to be tij. We prove that flat-bands 
emerge at E = —to. In this case the block matrix T in eq. (2) is written as 
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If we assume the eigenfunction <j)'j k with no amplitudes on B sites for the eigenvalue —to, the 
equation (T + £o-0( a ifcj a 2fc 5 ■ • ' > a ^fe)* = is transformed to M independent equations: 

t A 1 + t 12 A 2 + ■■■ + t 1M A M = (9) 

t 21 A 1 +t A 2 + -- - + t 2M A M = (10) 



where 



tM\A\ + tj\{2A 2 + h to^M 



• 1 • 



0. 



(11) 



(12) 



By adding eq. (4), we obtain N + M — 1 independent equations for z unknowns, which leads 
to z — N — M + 1-fold degenerate flat-bands at E = —to. The examples with to = 0, which 
result in flat-bands at E = 0, are shown in Figs. 3(b) and 4 (a). It is interesting to note that 
in the 3D example (Fig. 4(a)), another flat-band appears at E = —It. This is a unique feature 
of the model where the hopping amplitudes Uj between the subsets Si and Sj are independent 
of i and j, and M = 3. 

Finally, let us prove generalization (iv). Here, we maintain the original sites on the A 
sublattice, which has been eliminated on the partial line graphs. Even in this lattice, we may 
assume the wave function 0^, which has no amplitudes on B sites for the eigenvalue —t. 
Then we have N independent equations for z + \ unknowns, and can determine the z — N + 1 
independent solutions. Thus, there are flat-bands at E = — t with z — N + 1-fold degeneracy. 

In this paper, we have introduced a systematic method for constructing tight-binding 
models on generalized partial line graphs with flat energy bands. We have shown above only 
examples for the lattices divided originally into two sublattices. However, we can apply the 
method to the lattices with any number of sublattices. One realizes that the combinations 
of the described generalizations are numerous and useful for generating flat-band models. 
Since the flatness of the dispersion may lead to instabilities in the presence of interactions, it 
is important to study the many-body effects of the electrons on the generalized partial line 
graphs. For example, the ground state with a half-filled flat-band can be a partially polarized 
ferromagnetic state. It might be obvious that a ferrimagnetic state is a ground state in some 
of the generalized partial line graphs of Fig. 3(c), since one of them is the Lieb lattice, where 
it is proven rigorously that the ground state is a ferrimagnetic state. On the other hand, it 
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is well known that the fully polarized state on the flat-band is the unique ground state, 
when the flat-band is the lowest (highest) energy band and the single-particle density matrix 

Pxy ^ 

T c ytl^) * s irreducible. 3 Although the energy eigenvalue of the flat-band in the 
partial line graphs is not the lowest (highest), such a partially polarized ferromagnetic state 
can be a ground state when the band gaps on both sides of the flat-band exist for U/t <C l. 3 In 
fact, in a chain of square models, where a flat-band can be isolated, a flat-band ferromagnetism 
is realized for small U/t with appropriate parameters. 14 Since we can isolate a flat-band by 
tuning the parameters as shown in Fig2. 2(b) and 2(d), a partially polarized ferromagnetic 
state can be a ground state in certain generalized partial line graphs. 

We expect that the method is also useful in providing a guideline for synthesizing materials 
with flat-bands, since there is a large room for modifying the models. We can vary on-site 
energies, hopping amplitudes, and other parameters, while maintaining the flat-bands. Such 
flexibility is important in the view point of material designing. The fact that the on-site 
energies on B sites do not affect the flatness of the dispersion is particularly important. This 
fact indicates that the atomic species on B sites can be different from those in A clusters in a 
realistic system. Thus, we may expect that a generalized partial line graph may be realized in 
certain compounds. As a candidate, we expect that the quasi-one-dimensional lattice shown 
in Fig. 3(a), the 2D one in Fig. 3(b), and the 3D one in Fig. 4(a) will be realized in real 
materials. It is also useful that the position of the flat-band can be tuned by changing the 
hopping amplitudes, which indicates that the position of the flat-band can be moved toward 
the Fermi level without changing the carrier density by applying (chemical) pressure to the 
material. A clusters can be considered as z-fold degenerate orbitals as treated in refs. 11 and 
12, which is also useful in the view point of material designing. 

This work is supported by a Grant-in Aid from the Ministry of Education, Culture, Sports, 
Science and Technology of Japan. 
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